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Abstract 

The new method of solving quantum mechanical problems is pro- 
posed. The finite, i.e. cut off, Hilbert space is algebraically imple- 
mented in the computer code with states represented by lists of vari- 
able length. Complete numerical solution of a given system is then 
automatically obtained. The technique is applied to Wess-Zumino 
quantum mechanics and D=2 and D=4 supersymmetric Yang-Mills 
quantum mechanics with SU(2) gauge group. Convergence with in- 
creasing cut-off was observed in many cases well within the reach of 
present machines. Many old results were confirmed and some new 
ones, especially for the D=4 system, are derived. Extension to D=10 
is possible but computationally demanding for higher gauge groups. 
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1 Introduction 



Since the original conjecture of Banks, Fischler, Shenker and Susskind of 
the equivalence between M-theory and the D=10 supersymmetric Yang-Mills 
quantum mechanics (SYMQM) a lot of effort has been put to understand, 
and ultimately solve, the latter P|-|[T7[|. 

The general, not necessarily gauge, supersymmetric quantum mechanical 
systems have much longer history p, |^. Claudson and Halpern considered 
for the first time the gauge systems also well before the BFSS hypothesis 
In particular a complete solution of the D=2, N=2 SYMQM was given 
there (see also Ref. [Q) . Later the specific gauge models were studied in 
more detail and some candidates for the ground state were constructed in 
the, gauge invariant, Born-Oppenheimer approximation 0. 

Another important development was achieved by de Wit, Liischer and 
Nicolai who have shown that the spectrum of supersymmetric Yang-Mills 
quantum mechanics is continuous ||, ^ due to the cancellations between 
fermions and bosons. This is different from the pure Yang-Mills case where 
the transverse fluctuations across the vacuum valleys do not cancel and ef- 
fectively block the valley, resulting in the discrete spectrum of the 0-volume 
glueballs ||TU[. The continuous spectrum was first regarded as a setback, 
however it has turned out into a virtue with the advent of the BFSS hypoth- 
esis and new interpretation in terms of the scattering states. Still this new 
connection requires existence of the localized state at the threshold of the 
spectrum - the supergraviton This question triggered intense studies 

of the Witten index of SYMQM for various D and different gauge groups 



12| , 0, |14|, [1^. Powerful techniques were developed to calculate analytically 
nonabelian integrals [|TB], |I7|, related to the Witten index. They were ac- 
companied by complementary and original numerical methods |]19|, The 
emerging picture is rather satisfactory indeed: for D < 10 there is no thresh- 
old bound state while for D = 10 there is one, exactly as required by BFSS. 
This has been proven for N=2, but the evidence is being accumulated that 
it holds for higher N, as well as for other gauge groups. 

The large N limit of SYMQM, which is relevant to M-theory, was studied 



in the framework of the mean field approximation in Refs |21, p2l. This pro- 



vided an interesting realization of the black hole thermodynamics predicted 
by the M-theory [|Il], 



In spite of all above results, SYMQM remains unsolved for D > 2. There- 
fore it seems natural to study this model with lattice methods which proved 
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so successful in treating more complex field theoretical systems. Such a 
programme has been proposed in Ref. [Q, beginning with the yet simpler 



(quenched, D=4, N=2) case which was later extended to higher gauge groups 
2 < iV < 9 p5|]. The asymptotic behavior in N was observed, and indications 
of an onset of the interesting phase structure was found in agreement with 



Refs. pl|, p2| . Including dynamical fermions is possible for D=4, and may be 
feasible at D=10 for the first few N. However an arbitrary N case is plagued 
by the sign problem caused by the generally complex fermionic determinant. 

Another interesting approach studied by many authors, follows the Eguchi- 
Kawai trick to trade entirely the D-dimensional configuration space into a 



group space |£6|, Proceeding in this way one is led to consider the fully 
reduced (to a single point in the Euclidean space) model which nevertheless 
possesses a version of the superymmetry pSl . Many analytical and numerical 



results were obtained in this way (for a review see e.g. |^) . Numerical sim 



ulations were pushed to quite high N in simplified models However the 
sign problem which is also present there limits the Monte Carlo approach. 
This may be alleviated by the new method proposed recently in Ref. . 

It is important to remember that the subject has a lot of overlap with 
the small volume study of gauge theories where the valuable expertise has 



been accumulating for a long time |T^, Q |3J]. Although the final goal 



there is to increase the space volume and to match ultimately the standard 
large volume physics, the starting point is the pure Yang-Mills quantum 
mechanics identical to the bosonic part of our supersymmetric systems. In 
fact the classical results of Liischer and Miinster are the special case of the 
solution of the D=4 SYMQM in the gluino-free sector of the latter. This 
will be shown in Sec. 5. Recently van Baal has used the full machinery of 
the small volume approach to study the supersymmetric vacuum state of the 
more complicated, compact version of D=4 supersymmetric YM quantum 
mechanics . 

Summarizing, even though above quantum mechanics is much simpler 
that the M-theory, it remains unsolved and thereby still poses an interesting 
challenge. 

In this paper we propose a new approach and present a series of quanti- 
tative results for simpler systems not always solved up to now. To this end 
we use the standard hamiltonian formulation of quantum mechanics in the 
continuum^], construct explicitly the (finite) basis of physical states and cal- 

^This work was inspired by the discussion with C. M. Bender and the methods he 
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culate algebraically matrix representations of all relevant observables. This 
done, we proceed to calculate numerically the complete spectrum, the en- 
ergy eigenstates, and identify (super) symmetry multiplets. This approach 
is entirely insensitive to the sign problem and is equally well applicable to 
systems with and without fermions. Similarly to the lattice approach, the 
method has an intrinsic cut-off: any quantitative results can be obtained only 
within the finite dimensional subspace of the whole Hilbert space. It turns 
out, however, that in all cases studied below many important characteristics 
(in particular the low energy spectrum) can be reliably obtained before the 
number of basis vectors grows out of control. 

In the next section we describe the method and its algebraic computer 
implementation. In sections 3 and 4 the spectra of Wess-Zumino quantum 
mechanics and D=2 supersymmetric Yang-Mills SU(2) quantum mechanics 
are derived. Section 5 contains main results of this paper. It is devoted to 
D=4 supersymmetric Yang-Mills quantum mechanics with the SU(2) gauge 
group. The global picture of the whole system, with the quantitative spectra 
in all fermionic channels and their supersymmetric interrelations, will be pre- 
sented for the first time. Summary and discussion of the future applications 
follow in the last two Sections. 



2 Quantum mechanics in a PC 

We begin with the simple observation that the action of any quantum me- 
chanical observable can be efficiently implemented (e.g. in an algebraic pro- 
gram) if we use the discrete eigen basis 

{|n>}, |n >= ^(at)"|0 >, (1) 



of the occupation number operator a) a. For example, the bosonic coordinate 
and momentum operators can be written as 

x=^={a + a)), p= — i=(a — a'), (2) 

where all dimensionfull parameters are set to 1. Since typical quantum ob- 
servables are relatively simple functions of x and p, they can be represented 



developed in studying various quantum mechanical systems . 
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as the multiple actions of the basic creation and annihilation operator^. 
Including fermionic observables is straightforward and will be done in sub- 
sequent Sections. Generalization to more degrees of freedom is also evident 
and will be carried out separately for individual systems. 

The first step to quantify above considerations is to implement the Hilbert 
space in an algebraic program. Any quantum state is a superposition of 
arbitrary number, n^, of elementary states \n > 

\st >= S7^a/|n(^) >, (3) 

and will be represented Mathematica list f\ 

St = K, {ai, . . . , a„J, {n«}, {n^'^, {n^-)}}, (4) 

with Us + 2 elements. The first element specifies the number of elementary 
states entering the linear combination, Eq.(|^), the second is the sublist sup- 
plying all (in general complex) coefficients a/, / = 1, . . . , n^, and the remain- 
ing Us sublists specify the occupation numbers of elementary, basis states. 
In particular, convention implies for an elementary state 

|n>^{l,{l},M}. (5) 

In the case of more degrees of freedom the single occupation numbers n^^"^ 
in Eq.(^ will become lists themselves containing all occupation numbers 
corresponding to the individual degrees of freedom. 

In the next step we implement basic operations defined in the Hilbert 
space: addition of two states, multiplication by a number and the scalar 
product. All these can be simply programmed as definite operations on 
Mathematica lists. Table 1 displays corresponding routines which transform 
lists in accord with the principles of quantum mechanics. This done, it is 
a simple matter to define the creation and annihilation operators which act 
as a list-valued functions on above lists. This also defines the action of the 
position and momentum operators according to Eq.(|^). Then we proceed to 
define any quantum observables of interest: hamiltonian, angular momen- 
tum, generators of gauge transformations and supersymmetry generators, to 
name only the most important examples. 

^Thc method can be also extended to non polynomial potentials. 

•^Of course any other algebraic language can be used. In more advanced and time 
consuming applications we shall be using much faster, compiler based, languages anyway. 
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iiiapi)iiig 



any state 



\st > 

Sti > +\st2 > 

a\sti > 

< Sti\st2 > 



list 



number multiply 
scalar product 
empty state 
null state 



sum 



> 





{listi, list2) —>■ lists 
listi list2 
{listi, list2) — >■ number 



{0,{}} 



Table 1: Quantum states in the Hilbert space and their computer imple- 
mentation. 

Now our strategy should be obvious: given a particular system, we define 
the list corresponding to the empty state, then we generate a finite basis 
of Ncut vectors and calculate matrix representations of the hamiltonian and 
other quantum operators using above rules. Thereby the problem is reduced 
to a simple question in linear algebra, namely the spectrum of the system is 
given by the eigenvalues and eigenvectors of the hamiltonian matrix, their 
behaviour under rotations by the angular momentum matrix, etc. 

Of course the most important question is how much the ultimate physics 
at Ncut = oo is distorted by the finite cut-off. This can be answered quan- 
titatively by inspecting the dependence of our results on N^ut for practically 
available sizes of the bases. In all systems studied so far (and discussed be- 
low) the answer is positive: one can extract the meaningful (i.e. N^ut = oo) 
results before the size of the basis becomes unmanageable. However this can 
be answered only a posteriori and individually for every system. 

The method outlined above turns out to be a rather powerful tool capable 
to solve quantitatively various quantum mechanical problems with finite but 
large number of degrees of freedom. 

3 Wess-Zumino quantum mechanics 

This system is obtained by dimensional reduction of the four-dimensional 
Wess-Zumino model leaving only the time dependence of the dynamical de- 
grees of freedom. Resulting quantum mechanics (WZQM) has one com- 
plex bosonic variable (j){t) — x{t) -\- iy{t) = xi{t) -\- ix2{t) and two complex 
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Grassmann- valued fermions ipa{t), a = 1,2 The hamiltonian reads 



H = lipl+pD + \F\' + {{m + 2^?0)^iV'2 + h.c.) , (6) 



where the auxihary field F is determined in terms of by the classical equa- 
tion of motion, 

F = -{m<p + g<f)'), (7) 

since it enters only quadratically into the action. Independent variables sat- 
isfy the standard commutation relations 

[Xi,Pk] = iSik, {lpa,'ipl} = Sal3, (8) 

which allow us to parameterize these coordinates in terms of the creation 
and annihilation operators 



x = —{acc + al), = —^{a^ - al), (9) 



^a = /a, ^i = fl (11) 

which in turn satisfy 

4] = ^ik, {fa, 4} = Saf3. (12) 

Now we implement this quantum mechanical system in the computer as 
explained in the previous Section. The empty state will be represented by a 
list 

1(0, 0), (0, 0) >^ {1, {1}, {{0, 0}, {0, 0}}, (13) 

where the first brackets specify occupation numbers of the bosonic, and the 
second fermionic, states. To define properly the fermionic creation/annihilation 
operators we follow the original construction of Jordan and Wigner 



U = Ilp<a{-lf'cT- /t =n^<„(-l)^^a+, (14) 

where the fermionic number operator = f^fa (no sum) and are rising 
and lowering operators, commuting for different a, with (a^)^ = . The 
Jordan- Wigner phases ensure uniform anticommutation rules for fermionic 
operators. 
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We proceed to construct the finite eigen basis of tlie fermionic, Fa, and 
bosonic, Bi, number operators. For more tlian one degree of freedom the 
organization of a basis and definition of the cut-off, Ncut, is not unique. In this 
case we have chosen the cut-off to be the maximal number of bosonic quanta. 
That is, the basis subject to the cut-off Ncut consists of all orthonormal 
elementary states with B = B1 + B2 < Ncut and all allowed fermionic quanta, 
i.e. F = Fi + F2 < 2. In practice the basis is created by action on the 
empty state with all elementary independent monomials of bosonic creation 
operators up to Ncut th order, followed by the three independent monomials 
of fermionic creation operators. Hence the Hilbert space, subject to the 
cut-off Ncut, has 2{Ncut + l)(^cMt + 2) dimensions. 

Given the basis, the matrix representation of the hamiltonian (and any 
other observable of interest) can be readily calculated with our computer 
based "quantum algebra". The spectrum and energy eigenstates are then 
obtained by numerical diagonalization of the hamiltonian matrix. 

Figure ^ shows the spectrum of low energy states as a function of the cut- 
off up to Ncut=^5 for m = g = 1. Energies of the fourth and fifth level are 
shifted by 1 to avoid confusion of levels for low Ncut- Clear convergence with 
Ncut is seen, well within a capacity of a medium class PC. The convergence is 
faster for lower states. This general feature of our method can be easily un- 
derstood. The basis generated by creation operators, Eq.(|I2|), is nothing but 
the eigen basis of some normalized harmonic oscillator. It is obvious that the 
ground state can be easier approximated by a series of harmonic oscillator 
wave functions that the higher states with all their detailed structures, ze- 
roes, etc. Quantitatively, at A^c«t=10 the exact supersymmetric ground state 
energy is reproduced with the 3% accuracy which further improves to below 
1% at Ncut=^5 0- Supersymmetric pattern of the spectrum is also evident. 
First, the lowest bosonic state has no fermionic counterpart, its energy tends 
to zero, hence we identify it with the supersymmetric vacuum of the model 
Q. Second, all higher states group into bosonic-fermionic multiplets with the 
same energy. Splittings inside the multiplets decrease with Ncut as shown in 
Table 2. We conclude that supersymmetry at low energies is restored at the 
level of few percent in the cut Hilbert space with states containing up to 15 
bosonic quanta. 

^Since the exact value is zero, we take the first excited energy as a reference scale. 
^The model exhibits additional two-fold degeneracy in bosonic and fermionic sectors 
which is not related to supersymmetry |38| 
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Figure 1: Spectrum of WZ quantum mechanics as a function of cut-off, 



cut- 



Let us see how the above convergence to the supersymmetric spectrum 
shows up in the Witten index, which in this approach can be calculated 
directly from the definition 



Iw{T) = J:^{-lf'expi-TE,) 



(15) 



Fig.|^ shows Witten index for 4 < Ncut < 15. Nice convergence to the 
exact, time independent, result IwiT) = 2 is observed. Two things are 
happening which ensure this behaviour. First, there are two (due to the above 
degeneracy within bosonic and fermionic sectors) supersymmetric ground 
states. At the largest cut-off, Ncut=^5, their energies are not exactly zero 
but tiny Eq^ = 0.0046, Eq^i = 0.0071, hence giving rise to the exponential 
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Ncut 


Ml 


M2 


M3 


9 


0.1321 


0.1574 


0.1835 


10 


0.0924 


0.1213 


0.1245 


11 


0.0595 


0.0941 


0.0886 


12 


0.0420 


0.0657 


0.0609 


13 


0.0282 


0.0440 


0.0407 


14 


0.0212 


0.0332 


0.0287 


15 


0.0149 


0.0207 


0.0192 



Table 2: Splittings within the first three supersymmetric multiplets for dif- 
ferent cut-ofi^s. 

fall-off at large T with a very small slope. Second, cancellations between 
non zero energy states within supersymmetric multiplets are not exact and 
still leave exponential terms, albeit with smaller and smaller (with increasing 
Ncut) coefficients. This gives the residual time dependence which is however 
much weaker than the exponential one. 

Another, rather general feature of Iw{T)., can be also observed here. Ex- 
actly at T = Witten index vanishes for every Ncut just because the bases 
generated with our procedure automatically satisfy a "global" supersymme- 
try requirement. Namely the number of bosonic and fcrmionic states in a 
basis is the same for every Ncut- Therefore we conclude that 

lim Iw{^,Ncut)^^. (16) 

On the other hand exact supersymmetry requires that this complete corre- 
spondence between fermionic and bosonic states is violated since the ground 
state at ii^ = does not have its fermionic counterpart while each non zero 
energy state has. This apparent paradox is reconciled by noting that exact su- 
persymmetry emerges only in the limit of the infinite Ncut- In another words 
the unbalanced fermionic state is pushed (with increasing Ncut) to higher 
and higher energies and eventually "vanishes from the spectrum''. However 
it leaves a visible effect - Witten index has a discontinuity at T = which 
can be expressed as the noncommutativity of the T — > and Ncut ~^ 00 
limits 

hm lim7H^(r,7V^t)^lim hm Iw{T,N^t). (17) 

■'Vcut— »oo 1 — »0 1 — »0 A^cut— »oo 
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Figure 2: Witten index for Wess-Zumino quantum mechanics. 



The right hand side of the above equation, known as the bulk contribution, 
may not be integer for the continuous spectrum. 

Above resuhs were obtained for m = g = 1, but the method works equally 
well for any other choices of parameters. The case m = has a rotational 
symmetry and can be solved by other methods I^TJ. To our knowledge no 



quantitative calculation of the spectra and wave functions for m ^ exists. 
As mentioned earlier this approach provides a complete quantum mechanics 
of a system. In particular, wave functions in the coordinate representation 
are simply given by the well defined linear combinations of the harmonic oscil- 
lator wave functions. More detailed discussion of the model, and comparison 
with other authors will be presented elsewhere. 



4 D=2, supersymmetric, SU(2) Yang-Mills 
quantum mechanics 

This system, although solved analytically 0, has two new features: gauge 
invariance and continuous spectrum. Therefore we have chosen it as a next 
exercise. Reducing from D = 2 to one time dimension one is left with 
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the three real, colored, bosonic variables Xa{t) and three complex, fermionic 
degrees of freedom iJait) also in the adjoint representation of SU(2) , a = 
1,2,3. 

The hamiltonian reads 0] 

H = ^PaPa + igeabcijixblpc, (18) 

where the quantum operators V','?/'^ satisfy 

[Xa.Pb] = iSab, {^a, V'fe} = ^ab- (19) 

Hence they can be written in terms of the creation and annihilation operators 
as before 

Xa = ^{aa + al), p^ = -^{aa-al), (20) 

i^a = fa, = fl (21) 

To implement fermionic creation and annihilation operators we again used 
the Jordan- Wigner construction. 

The system has a local gauge invariance with the generators 

Ga = eabc{xbPc - ii'li^c)- (22) 

Therefore the physical Hilbert space consists only of the gauge invariant 
states. This can be easily accommodated in our scheme noting that the gauge 
generators of the SU(2) are just the angular momentum operators acting in 
color space. In fact the fermionic part of Eq.(p2D can be also interpreted 
in this way since the momentum canonically conjugate to ip: vr^ = iip^ . 
Therefore we construct all possible invariant under SU(2) combinations of the 
creation operators (referred for short as creators) and use them to generate a 
complete gauge invariant basis of states. There are four lower order creators: 

aWa = (aa), (23) 

4/i = («/), 

^abcaifbfl = (aff), 

eabcflfhl = (///)• (24) 
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Pauli principle implies that (//) = (a/)^ = {affY = {f ffY = 0, therefore 
the whole basis can be conveniently organized into the four towers of states, 
each tower beginning with one of the following states 

|0f>=|0>; |li.>=(a/)|0>, |2f >=(a//)|0>, |3^ >= (///)|0 >, 

(25) 

where we have labeled the states by the gauge invariant fermionic number 
F = flfa- The empty state (and its Mathematica representation) reads 

1(0, 0, 0), (0, 0, 0) >^ {1, {1}, {{0, 0, 0}, {0, 0, 0}}}, (26) 

with the obvious assignment of bosonic and fermionic occupation numbers. 
To obtain the whole basis it is now sufficient to repeatedly act on the four vec- 
tors. Eg . (|25|) , with the bosonic creator (aa), since action with other creators 
either gives zero, due to the Pauli blocking, or produces hnearly dependent 
state from another tower. Basis with the cut-off Ncut is then obtained by ap- 
plying (aa) up to Ncut times to each of the four "ground" states. Of course 
our cut-off is gauge invariant, since it is defined in terms of the gauge invari- 
ant creators. Moreover, since the gauge generators can be implemented into 
the algebraic program as any other observables, we can directly verify gauge 
invariance of the basis and any other state in question. 

Once the above basis is constructed, we proceed to calculate the matrix 
representation of the hamiltonian, the spectrum and the eigenstates accord- 
ing to Sec. 2 |. Again all computation can be done on a small PC and the 
longest run, corresponding to 20 bosonic quanta, takes approximately 500 
sec. It is worth to mention that the most time consuming are the algebraic 
operations on the states (c.f. Table 1) in the PC-based abstract Hilbert space. 
Consequently the program spends most of the time calculating the matrix 
representations of various operators [] while the numerical diagonalization is 
rather fast. 

The hamiltonian, Eq.(|I^) can be rewritten 

H = ^PaPa + gXaGa, (27) 

hence it reduces to that of a free particle in the physical, gauge invariant 
basis. It follows that it preserves the gauge invariant fermionic number and, 

^AU numerical results are tor g — 1. 

^In fact this part can be (and often is) done analytically in the integer arithmetic. 
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consequently, can be diagonalized independently in each sector spanned by 
the four towers in Eg. (|25|) . 

The spectrum is doubly degenerate because of the particle-hole symmetry 
which relates empty and filled fermionic states {\0f l-ip >) and their 1- 
particle 1-hole counterparts {\1f >^ |2i? >). This symmetry is not violated 
by the cut-off and indeed we see it exactly in the spectrum. On the other 
hand, supersymmetry connects sectors which differ by 1 in the fermionic 
number e.g. ( 10^ ><->■ \1f > ) and in general is restored only at infinite 
Ncut ■ A good measure of the SUSY violation is provided by the energy of 
the ground state which should be 0. We see that the energy of the lowest 
(doubly degenerate) state converges to but rather slowly, c.f. Fig.|^. This is 
interpreted as the indication that the spectrum is continuous at infinite cut- 
off. Indeed it is hard to approximate the non localized state by the harmonic 
oscillator states (this in fact is our basis) and consequently the convergence 
of such an approximation must be slow |. This deficiency can be sometimes 
turned into an advantage as will be seen in the D=4 case. 



0.25 



0.2 
Eo 0.15 
. 1 




0.05; ; 

6 8 10 12 14 16 18 20 

Ncut 

Figure 3: D=2 supersymmetric Yang-Mills quantum mechanics. The energy 
of the lowest state as a function of he cut-off and the 1/Ncut fit (dotted line). 

^This argument can be turned into a proof that the free spectrum converges hke 
1 / Ncut which is also confirmed here (c.f. Fig.^. 
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Surprisingly however, there exists a particular scheme of increasing the 
basis which renders exact supersymmetry at every finite cut-off N^ut in this 
model. To see that, let us first discuss tests of the SUSY on the operator 
level. SUSY generators read 

Q = i^aPa, Qa = '4'aPa, (28) 

with ijj = ip^ and satisfy the algebra 

{Q,Q} = 2H-gXaGa, = = 0. (29) 

Since the matrix elements of the SUSY generators can be easily calculated 
in our PC-based Hilbert space approach, one can readily check the matrix 
element version of Eqs.(P^. Indeed we find that these relations are satisfied 
to better and better precision with increasing N^ut- If so, it is natural to ask 
what is the spectrum of the hamiltonian matrix H^^ defined by Eq. ( [25| ) at 
finite cut-off. To this end we diagonalize the matrix 

H^%' = li< N\Q\M >< M\Q\N' > + < N\Q\M >< M\Q\N' >), (30) 

where A^, A^' and M label vectors of our finite basi^. Of course SUSY gener- 
ators mix the fermionic number, therefore one should combine all four towers 
of Eq.(PB[) into one big basis. One more trick is required to achieve exact 
SUSY at finite cut-off: we have to increase the basis allowing from the be- 
ginning for the disparity between the fermionic and bosonic sectors. That is, 
the size of the basis grows as: 2 + 4 + 4 + . . .. Hence dimension of the cut 
Hilbert space is 2 -|- ANcut in this scheme. With this choice the spectrum of 
the H^^ has exact superymmetry as shown in Table 3. 

As required, the ground state has zero energy and does not have the su- 
persymmetric image while higher states form supersymmetric doublets with 
the same energy. Additional degeneracy is caused by the particle-hole sym- 
metry as explained earlier. All non zero eigenvalues (with finite "principal" 
quantum number) tend to zero with increasing N^ut and form, at N^ut = oo, 
the continuum spectrum of a free hamiltonian, Eq. (|27|) . To reach the solu- 
tions with non zero energy, one should appropriately scale the index of an 
eigenvalue with Ncut ■ The four towers of eigenstates are in the direct corre- 
spondence with the four gauge invariant plane waves constructed in Ref.[0. 

^The second term in Eq.(G9h does not contribute in the physical basis. 
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Table 3: The spectrum and the degeneracy, g, of the hamiltonian H'^^ de- 
fined at finite Ncut as {Q, Q}/2. 

To our knowledge this is the first case where the cut-off preserves exactly 
the supcrsymmctry. On the other hand this success of the particular cut-off 
scheme may be related to the simple structure of the D=2 SYM quantum 
mechanics and its complete solubility. Nevertheless, it is important to keep 
in mind that at finite Ncut one is free to define the hamiltonian within the 
" 0(1 /Ncut)" tolerance, i.e. as long as various definitions converge to the 
same limit at infinite N^ut- Wc decided to use because it is the anti- 

commutator of charges which is used to derive basic properties of the SUSY 
spectrum. Further study of this exact realization of supersymmetry will be 
continued elsewhere. 

Finally let us shortly discuss the Witten index for this model. Because 
the particle-hole symmetry interchanges odd and even fermionic numbers, 
the Witten index vanishes identically. This is also true for any finite cut-off 
since Ncut preserves particle-hole symmetry. Nevertheless one can obtain a 
non trivial and interesting information by defining the index restricted to the 
one pair of fermion-boson sectors. For example 

7w^(r)(°'i) = E,,^,=o,i(-l)^^ exp i-TE,), (31) 

where the sum is now restricted only to the F = and F = 1 sectors. Since 
supersymmetry balances fermionic and bosonic states between these sectors 
(with the usual exception of the vacuum), the restricted index is a good 
and non trivial measure of the amount of the violation/restoration of SUSY, 
even when the total Witten index vanishes. Obviously 1^^'^^ = and 
^w^^ — ^w'^^ — ■ Studying restricted index is particularly interesting in 
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this model since, due to the continuum spectrum, it does not have to be 
integer and provides some information about the density of states. We have 
calculated the index from our spectrum of the original hamiltonian, Eq.(|T8]), 
for a range of cut-offs: 5 < Ncut < 20. Figure || shows the result for the 
2 + 4 + 4 + . . . scheme, i.e. when the basis is increased by every four vectors 
allowing from the beginning for the two unbalanced states from the empty 
and filled sectors each. We see a slow approach towards the time independent 
constant which seems to be 1/2. 



1 F 




^ 

1 2 3 4 5 

T 



Figure 4: Restricted Witten index for D=2 SYMQM. 

Witten index, as a sum over all states, provides an average measure of 
the breaking/restoration of SUSY. Indeed even though individual supersym- 
metric multiplets have not yet clearly formed (at cut-offs used for this calcu- 
lation), we see definite flattening of the T dependence of IwiT) which must 
then be a result of the average cancellation between many levels. Therefore 
it is easier to see the onset of the supersymmetric behavior in the Witten 
index than in the location of the individual levels. As before Fig. 4 shows 
that the limiting value of the index is discontinuous at T = with the same 
interpretation as in Sect. 3. However, contrary to the discrete spectrum of 
the Wess-Zumino model, the value of the index is not integer The value 

^"^It would be an interesting exercise to calculate the restricted index from the continuous 
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of the index at T = is the direct consequence of our scheme (2 + 4 + 4 + . . .) 
of increasing the basis and, of course, is scheme dependent. As a final re- 
mark we note an intriguing existence of a "universal" point in T at which 
the asymptotic value seems to be attained at all Ncut- 

5 D=4 supersymmetric SU(2) Yang-Mills quan- 
tum mechanics 

5.1 Preliminaries 

Dimensionally reduced (in space) system is described by nine bosonic coordi- 
nates x^(t), i = 1, 2, 3; a = 1, 2, 3 and six independent fermionic coordinates 
contained in the Majorana spinor « = 1? •••,4. Equivalently (in D=4) 

one could have chosen to impose the Weyl condition and worked with Weyl 
spinors. Hamiltonian reads 



H = Hb + Hp, (32) 

1 . ^ * J * J 

Hp = |e,,,^^^V6a;,^ (33) 

where ijj'^ is the transpose of the real Majorana spinor, and F in D=4 are 
just the standard Dirac a matrices. In all explicit calculations we are using 



Majorana representation of Ref.[43 



Even though the three dimensional space was reduced to a single point, 
the system still has the internal spin(3) rotational symmetry, inherited from 
the original theory, and generated by the angular momentum 

J" = e'^' - ^V^JS^Va) , (34) 

with 

S^'^ = -^[F,r^]. (35) 
Further, the system has gauge invariance with the generators 

Ga = take {xIpI - '-^Pl^?j , (36) 

spectrum of Ref . W . 
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and is invariant under the supersymmetry transformations generated by 



Qa = (rVa)aPa + t9eabc{^"'iJa)aXix',. (37) 

The bosonic potential (written now in the vector notation in the color space) 

V = ^S,fc(f^' X x'f, (38) 

exhibits the famous flat directions responsible for a rich structure of the spec- 
trum. At first sight one might expect that the spectrum of the purely bosonic 
(hence non supersymmetric) system does not have localized states because of 
these fiat valleys (c.f. Introduction). However the fiat directions are blocked 
by the energy of the transverse quantum fluctuations since valleys narrow as 
we move from the origin. As a consequence the spectrum of the model is 
discrete. Energies of the first lower states, known as zero volume glueballs, 
were first calculated by Liischer and Miinster p2]. On the other hand, in 



the supersymmetric system, transverse fluctuations cancel among bosons and 
fermions, valleys are not blocked, and the spectrum of the model is expected 
to be continuous 0]. The story has its interesting continuation in D=10 
model where the evidence of the threshold bound state was accumulated 



||13|| . Existence of such a state is necessary for the M-theory interpretation 
of the model where it is considered as a prototype of the graviton multiplet. 
Therefore detailed study of the low energy spectra of the whole family of 
Yang-Mills quantum mechanical systems, including identification of the lo- 
calized and non localized states is an important and fascinating subject. We 
shall face some of these questions also in the D=4 system. 

5.2 Creation and annihilation operators 

There are many ways to write quantum coordinates in terms of creation and 
annihilation operators 

K,at^] = 6:l, {f^,f^^} = 6^^, p,a = l,2. (39) 

The only constraint comes from the canonical (anti) commutation rela- 
tions 

[xlpt]=^5^'^5ab, {C,^f} = 5at. (40) 
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For bosonic variables we just use the standard extension of Eq.(PT|) to more 
degrees of freedom 

< = ^« + «?), v\ = ^{<-4)- (41) 

For fermionic variables we begin with the classical Majorana fermion in the 
Weyl representation P| 



^H/ = (C2,-cr,Ci,C2), (42) 

replace the classical Grassmann variables C? C* by fermionic creation and 
annihilation operators /, and transform to the Majorana representation. 
Final result is a quantum hermitean Majorana spinor 



/ 1 + ^ 



/ -fa-^f!+^f + f!J \ 
~^'^fa ~ fa ~ fa "I" ^/a 
~fa ~^ '^fa ~^ "^fa^ ~ fa 

V -^fl - fl + fV + ^/f ) 



(43) 



which satisfies Eq.(^0]) due to (^). Other choices of fermionic creation and 
annihilation operators are also possible |T5[. 

The next step is to define an empty state and its, computer based, alge- 
braic representation 

1(0, 0, 0), (0, 0, 0), (0, 0, 0), (0, 0, 0), (0, 0, 0) >^ (44) 
{1, {!}, {{0, 0, 0}, {0, 0, 0}, {0, 0, 0}, {0, 0, 0}, {0, 0, 0}}}, (45) 

where the first three vectors (in color) specify bosonic, and the last two 
fermionic, occupation numbers. 

The Jordan- Wigner transformation requires additional specification in 
this case. As before, the action of /* and on elementary states is defined 
by the spin-like raising and lowering operators corrected by the non local 
Jordan- Wigner phases 

fA = IiB<A{-lf''cr-B f\ = YiB<A{-lY^cj\. (46) 

However, since individual states are now labeled by a double index A, one 
must define the ordering of the two dimensional indices. We choose the 
lexicographic order. \i A = (a, a) and B = (b, j3) than 

(47) 

Any other unambiguous definition of the ordering is admissible. 
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5.3 Constructing the basis 

There are many ways to construct bases in higher dimensional systems. 
Apart from theoretical requirements one must also take into account practi- 
cal limitations of computer implementation. Of course the number of states 
at finite cut-off is generally bigger for higher D. However, the theory also 
has more symmetries and consequently the number of states in a particular 
channel can be kept manageable. 

First important symmetry is the conservation of the number of Majorana 
fermions 

[F,H]=0 F = r:^ra. (48) 

That is, the vector interaction Hp cannot produce Majorana pairs in this 
model. As a consequence the whole Hilbert space splits into seven sectors of 
fixed F = 0,1,..., 6. 

Second, the system has the particle-hole symmetry 

F (49) 

therefore the first four sectors F = 0, 1, 2, 3 contain all information about the 
spectrum. 

Further, the gauge invariance requires that our basis is built only from 
the gauge invariant creators. All this is not different from the earlier D=2 
case. 

The new element is the SO (3) invariance, with the angular momentum 
Eq.(0), which can be used to split further the problem into the different 
channels of fixed J. However, more we specify the basis more complicated 
its vectors become, and consequently more computer time goes into the gen- 
eration of such a basis and subsequent calculation of matrix elements. Sim- 
ilarly one might contemplate using the reduced version of Lorentz invariant 
composite operators as possible creators of a basis. Again this would gen- 
erate more complex basis since fields contain both creation and annihilation 
operators. Also, Lorentz covariance is not that relevant in our fixed frame, 
Hamiltonian formulation. 

Taking all above into account we have decided to produce the simplest 
basis of gauge invariant vectors using elementary creation operators. To this 
end we proceed as follows. Consider each fermionic sector separately, i.e. 
fix the fermionic number < F < 3. At given F create all independent 
vectors with fixed number of bosonic quanta B = a^^a^, and define a cut-off 
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as the maximal number of bosonic quanta B < Ncut, independently for each 
F. Hence Ncut can depend on F. To create all independent, gauge invariant 
states at fixed F and B consider all possible contractions of color indices in 
a creator of [F, B) order 

c^••<v.?^•c^ (50) 

for all values of the spatial indices i and a. All color contractions fall naturally 
into different gauge invariant classes. Creators from different classes differ 
by color contractions between bosonic and fermionic operators. For example 

(51) 

and 

(52) 

belong to different gauge invariant classes according to our definition. An- 
other example involves odd number of operators where one "contraction" 
consists of one triple of color indices coupled with the Levi-Civita symbol |^. 
Namely 

6e.e<^<ta^t„/t^n.t;;typt (53) 

and 

ec,ea^ai^at^a^^aff:^f,'^ (54) 

also belong to different gauge invariant classes. Different contractions within 
bosonic and fermionic family of operators do not lead to independent states 
hence are considered to be in the same class. Still there will be many lin- 
early dependent states generated by above creators with all values of spatial 
indices i and a. It is however a simple matter to recognize linearly inde- 
pendent ones given our rules of " quantum algebra" . We therefore leave it to 
the computer for the time being. At fixed F and B the final procedure is 
as follows: 1) identify all gauge invariant classes of creators, 2) loop over all 
values of spatial indices and for each zi, i^, cxi, cxj? create corresponding 
state from the empty state, Eq. (^5]) , 3) identify and reject linearly depen- 
dent vectors, 4) orthonormalize the remaining set of states. In principle this 
procedure depends exponentially on Ncut, hence can be further shortened 
and improved, however it is sufficiently fast for our present purposes. It also 

^^For the SU(2) gauge group, this step can be generahzed to higher SU(N) groups. 
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Table 4: Sizes of the bases generated in each fermionic sector, F. Ng is 
the number of basis vectors with given number of bosonic quanta, B, while 
E gives the cumulative size up to B. The last column gives the difference 
between the total number of the bosonic and fermionic states in all seven 
sectors. 

has an advantage of generating all channels of angular momentum (available 
within the cut-off Ncut)- This will prove useful in studying supersymmetry. 

Table IV shows the sizes of bases we have reached so far in each of the 
fermionic sectors. Due to the particle-hole symmetry the structure in the 
F = 4, 5, 6 sectors is identical with that in F = 2, 1, respectively. 

Calculation of the spectrum of H proceeds exactly as earlier. Naturally 
we can also readily obtain the spectrum of other observables, e.g. of the 
angular momentum, Eq.(0), as well as the angular momentum content of 
the energy eigenstates. For example the last row of Table 4 gives the highest 
angular momentum which can be constructed in each sector with the bases 
generated so far. 

5.4 Results 

The spectrum of the theory is shown in Fig.^. A sample of states is labeled 
with their angular momenta. It is important that our cut-off preserves the 
SO (3) symmetry. Consequently, the basis described in previous section con- 
tains complete representations of the rotation group for any Ncut ■ Hence 
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the spectrum displayed in Fig. 5 has appropriate degeneracies for each value 
of the angular momentum quantum number j. 

To maintain some clarity of the figure wc have cut arbitrarily the upper 
part of the spectrum, which extends to about Emax ~ 35 with current Ncut- 
One expects that the individual higher states have a considerable dependence 
on Ncut- However they also carry some relevant information, e.g. for quantum 
averages. 
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Figure 5: Spectrum of D=4 SYMQM. 

Apart from relating the bases in corresponding fermionic sectors, particle- 
hole symmetry also implies equality of all corresponding energy levels. Since 

our cut-off respects this symmetry we indeed observe, in sectors with F — 
6, 5, 4, exact repetition of the eigenvalues from F 0, 1, 2 subspaces. There- 
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fore only first four sectors are displayed in Fig.^. Evidently the spectrum of 
D=4 theory is very rich and raises many interesting questions. We shall dis- 
cuss some of them, beginning with the relation to the already known results. 



5.4.1 Correspondence with pure Yang-Mills quantum mechanics 



It turns out that the fermionic part of the hamiltonian, Eg. ([33|) , annihilates 
any state in the F = sector 

Hp\Op >= 0. (55) 

Therefore the supersymmetric hamiltonian, Eq.(^) reduces in that sector 
to the hamiltonian of the pure Yang-Mills quantum mechanics, well known 
from the small volume approach |T^, ^ This offers us a chance to 



compare one special case of the present spectrum with the classical results 
of Ref.|]52|. Figure ^ shows such a comparison for the first three states as a 



function of Ncut- Nice convergence of our results towards those of Liischer and 
Miinster is seen. They used 120 states in each angular momentum channel 
and variationally optimized frequencies of the basis states. Our F = basis 
contains 238 states which span all angular momenta up to j = 8, and has 
only 12 states with j = 0. No variational adjustment of the frequencies 
was attempted. In view of this we consider above agreement as a rather 
satisfactory confirmation of the present approach. The assignment of the 
angular momentum, which is done automatically within our method, is also 
the same. In particular, the characteristic to the 0-volume, ordering of the 
tensor glueball and the first scalar excitation is reproduced. 

5.4.2 The cut-off dependence 

The dependence of the spectrum on Ncut is of course crucial, as it provides 
the quantitative criterion if the size of the basis is adequate. It is displayed 
in Fig. 7 for the first few lower levels in all fermionic sectors. It is clear that 
the cut-off dependence is different in various fermionic sectors. For F = 0,1 
(and F = 6, 5) the dependence is rather weak, suggesting that the energies of 
lower levels have already converged (within few percent) to their asymptotic 
values. On the other hand, in F = 2, 3, 4 sectors convergence is considerably 
slower. As explained earlier, one can use the rate of convergence with Ncut to 
distinguish between the continuous and discrete spectrum. It was proved that 



in the case of the free spectrum the convergence is 0(1/ Ncut) i-e. slow |42| 
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Figure 6: First three energy levels in the F = sector of D=4 SYMQM for 
different cut-offs. Soild lines show the 0- volume results of Ref . [5^ . 



For the localized bound states we conjecture that the rate of convergence 
with Ncut is sensitive to the large x asymptotic of the wave function. We 
have observed this in a simple anharmonic oscillator and other models. This 
regularity is also clearly confirmed in the Wess-Zumino and D=2 SYMQM 
models discussed in previous Sections. 

Taking this into account we claim that the low energy spectrum of D=4 
SYMQM is discrete in F = 0, 1, 5, 6 and continuous in the F = 2, 3, 4 sectors. 
This is an interesting quantification of the result of Ref. which was men- 
tioned earlier. Since the fermionic modes are crucial to provide continuous 
spectrum, it is natural that it does not show up in the sectors where they do 
not exist at all, or are largely freezed out by Pauli blocking. 

Second result, which is evident from Fig. 7, concerns the supersymmetric 
vacuum in this model. Assuming that indeed the eigen energies have ap- 
proximately converged in F = 0, 1 sectors (none of them to zero) it follows 
that the SUSY vacuum can not be in empty and filled sectors {F = 1,5 is 
also ruled out by the angular momentum) The obvious candidates are 
the lowest states in F = 2, 4 sectors with their energy consistent with zero at 

-"^^In early attempts, empty and filled sectors were considered as possible candidates for 
SUSY vacuum. 
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Figure 7: Cut-off dependence of the lower levels of the D=4 SYMQM in 
four independent fermionic channels. 



infinite Ncut , c.f. Fig. 7. It follows, together with the conjectured correlation 
between localizability and Ncut dependence, that the SUSY vacuum in this 
model is non normalizable. Further evidence is based on the structure of the 
supersymmetric multiplets, and will be presented in the next section. Very 
recently van Baal studied F = 2 sector in a more complex, non compact, 
supersymmetric model of the same family Although our results cannot 
be directly compared, due to the specific boundary conditions required in 
||35|| , he finds that the energy of the lowest state in F = 2 sector is indeed 
very close to zero. The deviations from exact zero are caused by the SUSY 
violating boundary conditions. 

It is important to note that not all states in F = 2, 3, 4 sectors have to 
belong to the continuum. Supersymmetry together with the discrete spec- 
trum in F = 0, 1 channels implies existence of the normalizable states among 
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the continuum of F = 2,3,4 states. This will be discussed later, here we 
only give one explicit example. Indeed the energy of the F = 2, j = 1 state, 
shown in Fig. 7 (flatter curve beginning at E=6), has definitely weaker de- 
pendence on Ncut than the others. We interpret this clS db SI gnature that the 
lowest \2f, Ij > state is localized. This situation may be a precursor of the 
more complex phenomenon expected in the D=10 theory. There, the zero 
energy, localized bound state of Dq branes should exist at the threshold of the 
continuous spectrum. Present example suggests that one way to distinguish 
such a state from the continuum may be by the different Ncut dependence. 

5.4.3 Supersymmetry 

Operator level. To check the supersymmetry algebra 

{Qa, Q/3} = 26c.pH + 2gTipxlGa, (56) 

at finite cut-off, we sandwich both sides of this operator relation between 
our basis , introduce the intermediate complete set of states to saturate the 
product of Q's and test the equivalent relations between matrix elements 

< N\{Q^,Q(,}\N' >= 25,^ < N\H\N' > . (57) 

The gauge term in Eq.(|56D does not contribute since the basis is gauge in- 
variant at any Ncut- There is, however, one limitation. As is evident from 
Eq.(^)), supersymmetry generators change fermionic number by one, and 
can change number of bosonic quanta, B, at most by two 

<N\Qc.\M>^0 if Fm = Fm±1, \Bm-Bn\<2. (58) 

It follows that, if a cut off representation of unity \M >< M\, Bm < 
Ncut is used to saturate any matrix element of a product < N\Q\M >< 
M\Q\N' >, it will be biased for high states in N, N' sectors. Quantitatively, 
the super symmetric relations, Eq.(|57|) are valid if 

Bn , Bm> < Ncut - 2. (59) 

This limits the number of matrix elements available for the test. Still many 
interesting predictions for lower states can be verified and should be satisfied 
exactly at finite Ncut- 
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One more property of supersymmetry generators should be kept in mind. 
Since they change F by 1, they move between different fermionic sectors. 
Hence the generic matrix equation ([57|), when restricted to a particular sector, 
reads (a = no sum) 

< F\H\F' >= (60) 
< F\Q^\F -IxF- l\Qa\F' > + < F\Q^\F + 1><F + l\Qa\F' > . 

That is, the hamiltonian matrix receives contribution from a single fermionic 
sector, only if it is computed in empty or filled sectors. Of course above 
SUSY relations have many other consequences. Some of them will be used 
later, when discussing emergence of the supersymmetric multiplets. 

Within the above limitations we have checked some consequences of SUSY 
algebra at finite Ncut- The spectrum of the matrix agrees exactly with 
that shown in Fig.|^ and derived directly form the hamiltonian. Moreover, the 
off-diagonal commutators indeed vanish in our basis. Many other interesting 
questions emerge and will be studied elsewhere. 

SUSY on the level of states. The next goal is to identify supersymmet- 
ric multiplets in the spectrum and to asses the effect of the cut-off on their 
splittings and mixings. It is clear from Fig. 5 that any direct search for ap- 
proximately degenerate states is difficult and may not be conclusive. However 
our construction provides another simple way to move within the SUSY mul- 
tiplets. Namely it is sufficient to use the explicit matrix representation of 
supersymmetry generators in bases summarized in Table ^ If supersymme- 
try was exact, acting with Q on an eigenstate of H, would generate states 
within the same multiplet. At finite Ncut this is no longer true, however we 
expect that the resulting state, Ql'if > say , will be spread around the ap- 
propriate multiplet member (or members) in another fermionic sector. This 
information will then be correlated with the spectrum of Fig. 5. Important 
simplification comes from the rotational symmetry, which is exact at every 
Ncut in this approach. Supersymmetric generators carry the angular momen- 
tum 1/2 which clearly limits the range of possible targets (5|\E' >. Indeed, 
we confirm that 

<F,j|g|F±l,/>=0 if \3-3'\^\. (61) 

To proceed, we denote the fc-th energy eigenvector in the m-th fermionic 
sector and the ra-th angular momentum channel by \mF,nj;k >. The sub- 
script j will be omitted where evident. A summary of various transitions we 
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\\Qi\ 






1 > 1 


2 



have analyzed is shown in Fig.^ For example, the action of Qi on the first 
state in F = sector gives 

] \|2_lI^i„ i/o-oinJn„ n.-i \|2^ 

95%, 
(62) 

which provides rather satisfactory confirmation of the first-glance guess from 
Fig. 5. For higher states situation is more complex. The second scalar state in 
this sector goes under Qi into Qi\Of, 0^; 2 > which decomposes in 40% into 
the I li?, 1/2; 3,4 > doublet and 43% of yet higher, but also close in energy, 
l/2;5,6 > doublet. This is not surprising, since higher states have not 
yet converged to their infinite Ncut limit. 

Even more interesting is to analyze SUSY relations between F = 1 and 
F = 2 sectors. Again the lowest |li7,l/2;l > state transforms predomi- 
nantly (84%) into the third scalar \2f, Oj; 3 > which is closest in the energy. 
Another consequence of the SUSY algebra can be seen here. Naively one 
would expect that combining above transitions |0i?,0j;l >—>■ 1/2;1 > 
and \1f, 1/2; 1 >^ |2^,0j;3 > would lead to the connection of F = and 
F = 2 sectors and extension of the multiplet. This contradicts Eq. (|57|) which 
implies that the Ql is diagonal in the eigen basis of the hamiltonian. Solution 
of this apparent paradox is simple. There are two states with j = 1/2 and 
while the multiplet beginning in F = sector ends on one linear combina- 
tion, the second multiplet must start on the orthogonal one. In other words 
the supersymmetric images of |0i?,0j;l > and \2F,0j;3 > in F = 1 sector 
must be orthogonal. We confirm this exactly (e.g. to the full precision of 
Mathematica numerics)[^ 

< 2f,0,\Q\Qi\0f,0j >~ 10-'^ (63) 

The next (in the energy) |li?, 3/2; 1 > state goes under Qi in 30% to \2f, Ij] 1, 2, 3 > 
triplet, with the rest shared among higher F = 2 states. In fact, the 
|2ir, Ij-; 1, 2, 3 > states are those with the fast N^ut convergence, which we 
have already pointed earlier as good candidates for locahzed states. Very 
little of |li?,3/2; 1 > goes back to the F = sector (c.f. fat head arrows in 
Fig-D 

6.3%, (64) 



WQil 


1f,3/2;1 > II 


2 

F=0 


WQi 


|1f,3/2;1> 


2 



13 



This is because we use the limited basis Eq.(M), which guarantees Eq.(|57|) exactly. 
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since there is no state with similar energy there This means that the 
j = 3/2 states start a new SUSY multiplet which extends into F = 2 sector 
and not into F = one. 

We are using only one of the four SUSY generators just as an example. 
Similarly we do not consider systematically all members of j 7^ multiplets. 
The full analysis, which in fact would be quite straightforward in the present 
setup, would involve complete discussion of the extended, = 4, SUSY 
algebra. However it is out of the scope of this article and will be done 
elsewhere. 

All above examples have a common feature. Clear SUSY regularities show 
up if the states have already converged, i.e. if their energy depends weakly on 
Ncut ■ This is not the case for the lowest states in F = 2, 3, 4 sectors which, 
we believe, form a continuum at infinite cut-off. Accordingly, and similarly 
to the D=2 SYM QM case, identification of multiplets and the SUSY vacuum 
in this sector is not clear at the moment. As one possible signature of the 
vacuum one might take the average norm of the supersymmetric image of 
the lowest state | |(5i|2j7, 0^; 1 > |p/(iim[3i?]. Indeed it is smaller than that 
for other states and falls with Ncut- However it is too early to draw definite 
conclusions. 

Clearly we have only begun. The aim here is just to introduce the new 
approach and its potential leaving specific applications for more focused 
articles. 

Witten index. The total number of bosonic and fermionic states is the 
same for each B, c.f. last column of Table 4. This is not a direct consequence 
of supersymmetry, but rather of a combinatorics of binary fermionic systems, 
which is not spoiled by gauge invariance. Total number of states gives the 
Witten index at T = 0. Hence we expect 

^lim IwiO,Ncut)=0. (65) 

As was already discussed Witten index is discontinuous at T = 0. In par- 
ticular /vk(0) depends on the regularization, therefore this number is tied to 
our B - F symmetric scheme of increasing the basis. 

Figure 9 shows the (euchdean) time dependence of the Witten index 
calculated from our spectrum with up to five bosonic quanta. Clearly we 

^*A11 surrounding states have already converged within few percent which suggests that 
the above small rate is also a 0{1/ Ncut) effect. 

^^We know of no other method capable to provide such a detailed and quantitative 
characteristics as discussed in this Section. 
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Figure 8: Supersymmeric images of a sample of eigenstates. 



are far from the convergence, but some interesting signatures can be already 
observed. First, one can see how the discontinuity at T = emerges even at 
this early stage. Second, the exponential fall-off common to all curves at large 
(T > 5) times is evident. This is where the lowest state dominates. Since at 
this Ncut its energy is not zero, the fall-off is exponential. Finally, and most 
interestingly, we observe the flattening shoulder appearing at T ~ 2 — 3. This 
is the signal of the supersymmetric cancellations which occur on the average, 
even though the exact SUSY correspondence between individual states does 
not yet appear. The behavior with N^ut is not inconsistent with the exact 
bulk value 1/4 [|13| obtained also from the nonabelian integrals [16 . 
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Figure 9: Witten index of the D=4, supersymmetric Yang-Mills quantum 
mechanics for the number of bosonic quanta bounded by Ncut = 2, 5. The 
bulk value 1/4 is also shown. 

Clearly higher Ncut are needed, and, what is important, they are perfectly 
within the range of present computers. 

6 Summary 

The new approach to quantum mechanical problems is proposed. The Hilbert 
space of quantum states is algebraically implemented in the computer code. 
In particular realization, used here, states are represented as Mathematica 
lists. All basic operations on quantum states are mirrored as definite opera- 
tions on above lists. Any quantum observable is represented as a well defined 
function on these lists. This allows for automatic calculation of matrix rep- 
resentation of a hamiltonian and other quantum operators of interest. To 
this end we use the discrete eigen basis of the operators of the numbers of 
quanta. The length of above lists is not fixed. Similarly to the length of 
an arbitrary quantum combination of basic states, it can vary dynamically 
allowing for any number of quanta. 

Of course in any finite computer the maximal length of a combination 
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must be limited. Therefore we impose a cut-off N^ut which bounds from 
above the number of allowed quanta. The stringent and quantitative test of 
this approach is provided by checking the dependence of any physical observ- 
able, e.g. the spectrum or the wave functions, on the cut-off. This is similar 
to studying the infinite volume statistical systems via the finite size scaling. 
We have applied above technique to the three progressively more compli- 
cated systems: Wess-Zumino quantum mechanics, supersymmetric Yang- 
Mills quantum mechanics in D=l-t-l dimensions, and D=4 SYM quantum 
mechanics, both for the SU(2) gauge group. In distinction from many other 
approaches (for example lattice simulations) the method is completely insen- 
sitive to the sign problem and works equally well for bosonic and fermionic 
systems. 

For Wess-Zumino quantum mechanics we can calculate the discrete spec- 
trum for any values of parameters. Clear restoration of the supersymmetry 
was observed for the cut-offs well within the capabilities of present comput- 
ers. Witten index was also obtained and its convergence to the known result 
IwiT) = 2 is clearly seen. 

The next system, D=2 SYM QM, possesses the gauge invariance which 
was readily incorporated in our approach. The physical subspace of gauge 
invariant states was explicitly constructed. Known structure of the solu- 
tions in terms of four fermionic sectors was reproduced. Convergence of the 
method, and emergence of the supersymmetry, was also studied in this more 
difficult case of the continuum spectrum. Witten index, restricted to one 
supersymmetric branch of the model, was defined and computed for the first 
time. Clear, albeit slow, convergence to the time independent fractional value 
/^(T) = ±1/2 was observed. Moreover, we have found a special scheme of 
increasing the basis such that the supersymmetry is exact at any finite cut-off. 
This however, may be related to the exact solubility of the model. 

Finally, the method was applied to the unsolved up to now D=4 SYMQM. 
This much richer system has the SO (3) rotational symmetry inherited from 
its space extended predecessor. Our approach preserves this symmetry ex- 
actly at any finite cut-off. The Hilbert space splits again into seven sectors 
with fixed fermionic number. We have obtained the complete spectra in all 
these sectors and studied their cut-off dependence. The spectrum in F = 
sector agrees with the classical 0-volume calculation for pure Yang-Mills 



quantum mechanics ||32[. An efficient method to distinguish between the 
discrete and continuous spectrum was proposed. It turns out that the asymp- 
totics of the wave function at large distances determines the convergence of 
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our calculations with the number of allowed quanta N^ut- The continuous 
spectrum with non localized wave functions converges slowly (~ 0{1/Ncut)i 
while discrete, localized bound states lead to faster (sometimes even exponen- 
tial) convergence. Accordingly, we have found an evidence that the spectrum 
of D=4 SYM QM is discrete in F = 0,1,5,6 sectors while it is continuous 
in the F = 2,3,4 sectors. This provides an explicit realization of the claim 
of Ref. PI in particular fermionic sectors. Interestingly, localized states exist 
also in the sectors with continuous spectrum. This is a simple consequence of 
supersymmetry whose transformations move between adjacent fermionic sec- 
tors. Our method allows to monitor directly the action of SUSY generators 
and analyse supersymmetric images of any state. In this way we have iden- 
tified some candidates for lowest supersymmetric multiplets. They do not 
have the same energy at current values of the cut-off, however the splittings 
are small and consistent with vanishing at infinite Ncut- In particular, the 0- 
volume glueballs found in F = sector are relevant to fully supersymmetric 
theory in that there exist gluino-gluon bound states with the same masses. 
We have also found a candidate for the supersymmetric vacuum which seems 
to belong to the continuum. However identification of SUSY multiplets in 
the continuum part of the spectrum requires higher Ncut- Supersymmetry 
was also tested on the operator level. In particular we confirmed that the 
spectrum of the hamiltonian coincides with that of Q^. As a final application 
we have calculated the Witten index for this theory. It still depends strongly 
on Ncut ■ Nevertheless an early evidence for some of its asymptotic properties 
can be already seen. 



7 Future prospects 

The main goal of this work is to asses the feasibility of attacking with the new 
approach the BFSS model of M-theory. Certainly the method is more efficient 
for smaller number of degrees of freedom. Application to the Wess-Zumino 
quantum mechanics and D=2 Yang-Mills quantum mechanics give quite sat- 
isfactory and quantitative results including the new intriguing scheme which 
preserves exact supersymmetry for finite cut-off. For more complex, and 
correspondingly richer, D=4 SYMQM the method is able to provide new 
results including detailed information about the supersymmetric structure 
of the spectrum and the observables. Obviously we also see a room for im- 
provement which is especially needed in the continuum sector. However, and 
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that we find most important, further increase of the size of the Fock space 
is possible within the available technology. The whole programme can be 
(and is being) implemented in the standard, compiler based languages which 
usually improves performance by a factor 10-100. Further, the action of the 
quantum operators can be optimized taking into account symmetries of the 
states. Finally, one can go for more powerful computers. 

Taking all above into account, one can reasonably expect substantial im- 
provement in the quality of the present D=4, N=2 results. At the same 
time one should be able to study D=4 systems with higher N. To answer 
the main question: yes, we think that the quantitative study of the D = 10 
theory is feasible, and reaching current quality results for the D=10 is realis- 
tic, beginning with the SU(2) gauge group. D = 10 hamiltonian and spin(9) 
generators do not conserve fermionic number 0, The reason for this 

complication should be better understood in the first place. 

Apart from increasing the number of degrees of freedom, higher gauge 
groups pose an interesting problem of constructing gauge invariant states. 
This has already been done in the small volume approach for SU(3) and 
should not present any fundamental problem for higher N as well. At the 
same time the possibility of some large N simplification, specifically within 
the present approach, should be investigated. 

Last, but not least, we would like to mention a host of applications of this 
method to other quantum mechanical systems. For example one may sim- 
ply extend the 0-volume calculations to the full non supersymmetric QCD 
with dynamical quarks in fundamental representation. Apart from known 
glueballs, this would give us a spectrum of quark- made-hadrons in the "fem- 
touniverse" . 

As another rather different application, we mention that this program is 
already being used as a routine in solving, to a high precision, quantum me- 
chanics of a simple two dimensional building blocks of a prototype quantum 
computer [^. In particular, a complete quantum evolution in time with the 
time dependent hamiltonian, was straightforward to simulate. 

The D=4, N=2 SYMQM studied here has 15 degrees of freedom. There 
are many unsolved quantum mechanical systems with this or smaller com- 
plexity. Present approach should be well applicable in some of these cases. 

Obviously there are many routes which can be followed from this point 
and we are looking forward to explore some of them. 
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